On the existence of supersolid "^He monolayer films 
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Extensive Monte Carlo simulations of ^He monolayer films adsorbed on weak substrates have been 
carried out, aimed at ascertaining the possible occurrence of a quasi- two-dimensional supersolid 
phase. Only crystalline films not registered with underlying substrates are considered. Numerical 
results yield strong evidence that ''He will not form a supersolid film on any substrate strong enough 
to stabilize a crystalline layer. On weaker substrates, continuous growth of a liquid film takes place. 



PACS numbers: 75.10.Jm, 05.30.Jp, 67.40.Kli, 74.25. Dw 



INTRODUCTION 



The search for a supersolid phase of ^He is still on- 
going, not only in the bulk crystal, but also in films of 
helium adsorbed on different substrates. Indeed, long be- 
fore the observation of non-classical rotational inertia in 
solid *He by Kim and Chan^i it was proposed by Crow- 
ell and Reppy^ that the second layer of a '^He film ad- 
sorbed on graphite might display simultaneously super- 
fluid and solid behaviour; this contention has been re- 
cently re-iterated<^ The most recent, accurate theoretical 
studies cast serious doubts on this claim;4.!§ furthermore, 
one should note that the denomination "supersolid" is 
scarcely appropriate for such a system, as the (alleged) 
crystalline order of the second layer would not arise spon- 
taneously as required by the definition, but rather be 
induced by an external agent, namely the graphite sub- 
strate. More generally, any claim of supersolid behaviour 
made in the context of a quantum film, requires that the 
crystalline phase of the latter be incommensurate with 
the substrate upon which the film is adsorbed. 

The idea that, in a many-body system wherein parti- 
cle motion is almost entirely confined to two dimensions, 
quantum fluctuations in the transverse direction may give 
rise to a superfluid response, even if the effective two- 
dimensional (2D) particle density corresponds to a crys- 
talline phase in the strictly 2D system, has already been 
extensively explored in the context of the search for a 
superfluid phase of parahydrogen. However, flrst princi- 
ples calculations carried out for various substrates^ii did 
not yield any evidence of superfluidity, a fact that could 
be certainly ascribed to the strong tendency of parahy- 
drogen to solidify, which renders it qualitatively different 
from helium. It is therefore not implausible that the more 
quantal helium may display a different behaviour in sim- 
ilar physical conditions. 

The possible supersolid behavior of quasi-two- 
dimensional helium, with motion of atoms in one direc- 
tion restricted by a confining potential, was investigated 
in a recent theoretical study, in which a simple harmonic 
model of transverse confinement was adopted.?^ The re- 
sults of that study were interpreted by the authors as 
suggestive of a novel phase, simultaneously displaying a 
superfluid response and broken translational invariance. 



underlain by large atomic excursions in the transverse 
direction. 

As the harmonic approximation may reasonably well 
describe the potential experienced by helium atoms in 
the topmost layer of a thin film, or in the first adsorbed 
layer on a weak substrate, the question naturally arises 
of whether there exists an actual experimental system 
that could give rise to observable quasi-two-dimensional 
supersolid ^He. Viable candidates are the second incom- 
mensurate solid ''He layer on graphite;^ or perhaps the 
first layer on a weak substrate, barely strong enough to 
stabilize it; the latter scenario is explored in detail here. 

This paper illustrates results of Quantum Monte Carlo 
(QMC) simulations of an ensemble of N ^He atoms, re- 
garded as point particles, moving in the presence of an 
infinite, smooth planar substrate of area A, positioned 
at z=0. The nominal ''He coverage is 6* = N /A. The 
quantum-mechanical many-body Hamiltonian is the fol- 
lowing: 
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Here, A=6.0596 KA^, V is the potential describing the 
interaction between two helium atoms, only depending 
on their relative distance, whereas U is the potential de- 
scribing the interaction of a helium atom with the sub- 
strate, also depending only on the distance of the atom 
from the substrate. We use the accepted Aziz potentialii 
to describe the interaction of two helium atoms. 

The interaction of a helium atom with the substrate 
(i.e., the U term in ([T])) is described here by means of 
the so-called "3-9" potential: 






(2) 



which is a functional form obtained by integrating the 
Lennard- Jones potential over a semi-infinite, continuous 
slab. Here, D is the characteristic depth of the attractive 
well of the potential, whereas a is essentially the distance 
from the plane of the minimum of such a well; the po- 
tential is strongly repulsive for z < a, weakly attractive 
for z > a. 



The aim of Us-g is that of providing a fairly reahs- 
tic description of the environment experienced by a he- 
Uum atom in the vicinity of a relatively weak substrate. 
Substrate corrugation, which is physically responsible for 
the appearance of low-coverage insulating films registered 
with the crystalline structure of the underlying substrate, 
is purposefully omitted in our model, as we are only in- 
terested in incommensurate solid films. 

For comparison purposes, we also performed calcu- 
lations using the same harmonic potential of Ref. |^, 
namely 
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defined in terms of the root-mean-square excursion a 
away from the equilibrium position (located at Zo above 
the plane) of the physisorption potential experienced by 
a free ^He atom near a given substrate. Obviously, a is 
the only parameter in this model, the actual value of Zo 
being immaterial. 

The low-temperature thermodynamics of the many- 
body system described by Eq. lU has been investi- 
gated by QMC simulations based on the continuous-space 
worm algorithm, a methodology that has proven quite 
effective in the study of Bose systems. Details of the 
computational methodology adopted in this work are de- 
scribed in the literatureJ^^— The results presented here 
are based on simulations of systems comprising iV=144 
atoms; the imaginary time step utilized in most of the 
calculations is r = 1.5625x10"'^ K~^, which has been 
empirically found to provide numerical estimates indis- 
tinguishable, within statistical errors, from those extrap- 
olated to the r = limit. The substrate geometry is 
chosen rectangular in all simulations, so as to accommo- 
date a triangular crystal. Periodic boundary conditions 
are assumed in all directions, but the size of the system 
in the z direction (i.e., perpendicular to the substrate) is 
chosen large enough to make boundary conditions unim- 
portant. This set-up is standard for QMC simulations of 
adsorbed films J^ 

Our main result is that no quasi-two-dimensional su- 
persolid *He phase arises in the physical conditions de- 
scribed above. Specifically, strong atomic transverse con- 
finement (corresponding either to a deep attractive well 
D in C/a-g or to a low value of a in Uh) leads to the 
formation of non-superfluid crystalline monolayer films, 
whereas on weaker substrates continuous growth of su- 
perfluid (liquid) film is observed. No intermediate ho- 
mogeneous phase, simultaneously displaying superfluid- 
ity and broken translational invariance, is ever observed 
in our simulations, including any metastable "glassy" su- 
perfluid phase. These results and conclusion are largely 
consistent with those of Ref. |8|; in our view, however, 
any interpretation of the available numerical evidence in 
terms of possible supersolidity, is baseless. 

In the remainder of this manuscript we illustrate our 
results, discussing first those for for the more realistic 
[/3_g, and then the ones yielded by the simple confining 



TABLE I: Values of the coefficients a and D of the "3-9" po- 
tential describing the interaction between a helium atom and 
four metal substrates, listed in order of decreasing attractive- 
ness (i.e., the well depth D). 



Substrate | Reference 


a(A)|D(K) 


Al 


16 


2.96 


60.3 


Ni 


16 


3.26 


48.7 


Mg 


17 


3.26 


35.6 


Li 


il 


3.76 


17.9 



potential Uh] as we shall see, there exist no substantial 
qualitative differences between the physical behaviours 
observed in two cases. 



II. 



'3-9" POTENTIAL 



Table |T] shows values of the coefficients a and D which 
have been proposed (and utilized in previous works) to 
describe the interaction of a helium atom with various 
metal substrates r^iii based on the form ^ of the ph- 
ysisorption potential. The most attractive of them is 
Al, for which the well depth D is 60.3 K, the least at- 
tractive is Li {D = 17.9 K). R should be stated at the 
outset that, aside from the very simple form of the "3- 
9" model potential, these coefficients are only known to 
a typical accuracy of 10-20%. The purpose of using a 
simple model such as ^, even with the current uncer- 
tainty in the values of the coefficients, is that of gaining 
fundamental understanding on the physics of adsorbed 
monolayer and guiding the experimental search toward 
systems that may display novel behaviour. 

Calculations making use of the "3-9" potential were all 
carried out at the same coverage 9 = 0.08 A~^, chosen 
high enough to fall unambiguously into the crystal re- 
gion of the two-dimensional (2D) phase diagram of ^He 
at low T,^* but also sufficiently low to make transverse 
zero-point motion significant. We have carried out calcu- 
lations at low T, the lowest being T = 0.1 K. In general, 
results obtained at T < 0.5 K can be regarded as essen- 
tially ground state estimates, as no significant tempera- 
ture dependence can be observed at lower T. 

Figure [1] shows computed '*IIe density profiles n{z) in 
the direction (z) perpendicular to the substrate, for the 
four metal substrates listed in Table HI Specifically, 



(4) 



"•(^) = (IM) / dx dy p{x,y,z) 



where p is the three-dimensional ^He density. 

The most weakly attractive of the four substrates, 
namely Li, is the weakest known substrate on which ^He 
is predicted to form a stable superfluid monolayer, at 
a coverage of approximately 0.056 A~^. On increasing 
coverage, promotion of ^He atoms to second layer is ob- 
served before monolayer crystallization ; ^^i^*^ as seen in 
Fig. [T] from the broadly extended shape of n{z), which 
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FIG. 1: (Color online). *He density profile n{z) (in A"^) of 
He in tlie direction perpendicular to the substrate, computed 
at r=0.2 K for the four substrates listed in Table U The 
nominal coverage is 0=0.08 A~^ in all cases. 
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FIG. 3: (Color online). One-body density matrix n{r), com- 
puted at T=0.2 K for Li (triangles) and Ni (circles) substrates. 
The nominal coverage is 61=0.08 A~^, whereas the tempera- 
ture is T=0.2 K. No significant change is seen at lower T. 
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FIG. 2; (Color online). Reduced ^He pair correlation func- 
tion, defined as in Eq. ([5}, at r=0.2 K for three of the four 
substrates listed in Table |T] (result for Al is indistinguishable 
from that for Ni, on the scale of the figure). The nominal 
coverage is 0=0.08 A~^ in all cases. Also shown (dotted line) 
is the pair correlation function for two-dimensional *He at the 
same temperature and coverage. 



features a shoulder beyond its main peak at z « 4.5 A. 
The same physics is observed on a Mg substrate, only 
slightly more attractive than Li. On Ni and Al, on the 
other hand, n{z) displays a single peak, suggesting that 
^He atoms form essentially a monolayer, positioned closer 
to the substrate for the more attractive Al, with signifi- 
cant zero-point motion of atoms in the z direction. 

The 2D character of ''He films on Ni and Al can be 
established by looking at the the angularly averaged, re- 
duced pair correlation function g{r), with r = \/ x^ -I- y^ 
and 



9^^^y) = -Jffl \ d'X' dy' n{. 



x + x\y + y') n{x',y') (5) 



with n{x,y) = J dz p{x,y,z). The more closely g{r) 
mimics the pair correlation function of a strictly 2D "^He 
system of the same coverage, the closer the film is to 
being essentially 2D. Figure [2] shows g{r) for a "^He film 
of coverage 0, adsorbed on Li, Mg and Ni substrates, 
compared to the pair correlation function of '^He in two 
dimensions. As mentioned above, in two dimensions *He 
at this coverage forms an insulating (i.e., non superfluid) 
triangular crystal at low T, and that is the physics ob- 
served on the two stronger substrates, namely Al and Ni, 
for which g{r) closely approaches the 2D shape. Atomic 
zero-point motion in the direction perpendicular to the 
substrate only causes a small depression of the main 
peak. On the other hand, on the weakly attractive Li 
and Mg, promotion of ''He atoms to the second layer is 
shown by the finite value that g{r) takes at the origin, 
as well as by the less pronounced first peak and oscilla- 
tions. On these two substrates, no 2D crystal of helium 
forms, and the ''He superfluid fraction is essentially 100% 
at T ^ 0.5 K, whereas no detectable superfluid signal is 
observed on Ni and Al, down to a temperature T=0.1 
K. These two very different behaviours are also reflected 
in the one-body density matrix n(r) (shown in Fig. [3|), 
which displays an exponential decay for the film adsorbed 
on Ni and Al, whereas on Mg and Li its expected much 
slower power-law decay is barely detectable within the 
precision of our calculationi^i 

Looking at Table HI we see that the value of a for Mg 
and Ni substrates is the same, the different physics of ad- 
sorbed helium films arising exclusively from the different 
well depths (i.e., D). The question can be posed, there- 
fore, of whether on a substrate characterized by that very 
same value of a and one of D between those of Mg and Li, 
a "supersolid" He monolayer may exist, namely an ef- 
fectively 2D crystalline film, featuring a finite superfluid 
response, possibly within a narrow range of temperature. 
If one could establish that supersolid behaviour can be 
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FIG. 4: (Color online). Integrated pair correlation function 
computed at T=1.0 (solid line) and 0.5 (dashed line) K, for a 
^He film of coverage 0=0.080 A~^ on a "3-9" substrate with 
D=40 K and a=3.26 A. 
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FIG. 5: (Color online). One-body density matrix n{r), com- 
puted at T=1.0 K (triangles), T=0.5 K (circles) and T^O.25 
K (boxes), for a "'He film of coverage 61=0.080 A"^ on a "3-9" 
substrate with D=40 K and a=3.26 A. Statistical errors are 
of the order of the sizes of the symbols. 



observed at an attainable low T a narrow range of values 
of D, then issue would of course be that of fashioning 
such a substrate in the laboratory. One possibility might 
be that of "weakening" an attractive substrate by pre- 
plating it with a rare gas, such as Kri^ 

We have addressed this issue by performing the same 
calculations for a "fictitious" substrate of well depth D 
in the range between 36 and 40 K and with a=3.26 A, in 
search of a "sweet spot" wherein quantum oscillations in 
the transverse direction may give rise to novel physics. As 
shown below, however, no other phase is observed other 
than a (non-superfluid) solid and a (superfluid) liquid 
film. 

Figure |4] shows the integrated pair correlation function 
g{r) computed on such a model substrate, with £'=40 
K, at two low temperatures, namely £=1 K and T=0.5 
K. It is worth mentioning that the result at £=0.75 K 
is indistinguishable, within the statistical errors of the 
calculation, from that at T=0.5 K. Clearly, the physi- 
cal behaviour of the system changes qualitatively as the 
temperature is lowered. At the lower T, the pair correla- 
tion function features the regular oscillations proper of a 
crystal, and is essentially 2D in character; at the higher 
temperature, on the other hand, little or no structure 
appears beyond the first peak, and the result is consis- 
tent with liquid-like behaviour. Moreover, in the latter 
case the small but finite value of g at the origin points to 
thermally activated promotion of helium atoms to second 
layer, i.e., the quasi-2D crystal melts into a superfluid at 
a temperature £„ « 1 K. 

The superfluid character of the fllm at higher temper- 
ature can be established by direct computation of the 
superfluid density, by means of the well-known "winding 
number" estimator, as well as through the behaviour of 
the integrated one-body density matrix, shown in Fig- 
ure \5\ While at the higher T the function displays a 
slow decay (and correspondingly the value of the super- 



fluid density is approximately 50% for the finite system 
simulated), at lower T it decays exponentially, with no 
appreciable dependence on T, as shown in Figure [3 In 
fact, the result at £=0.25 K is indistinguishable from 
that at £=0.5 K, within statistical errors. 

The above results represent the general trend observed 
on performing calculations with any value of D, for fixed 
a. Increasing (decreasing) D has simply the effect of rais- 
ing (lowering) the temperature £„ at which the quasi-2D 
crystal melts. If D is large, the crystal melts into a nor- 
mal liquid, whereas superfluidity is observed for interme- 
diate values of D. On the other hand, as D approaches 
a value close to ^ 36 K, i.e. near that corresponding 
to a Mg substrate, then £„ -^ 0, i.e., no crystal film 
is observed. Simultaneous occurrence of superfluid and 
crystalline properties is never observed. 

Thus, the physical conclusion of this theoretical exer- 
cise, is that quantum-mechanical particle fluctuations in 
the transverse direction do not lead to the stabilization 
of a quasi-2D supersolid phase of "^He, at least within the 
framework defined by a potential such as ^ . 



III. HARMONIC POTENTIAL 

The conclusion at which we arrived in the previous 
section would appear to clash with recent predictions of 
supersolid behavior from Cazorla et al. (Ref. ^), who 
performed a strictly £=0 calculation based on a har- 
monic model of confinement (Eq. [3|). It seems difficult 
to ascribe such a different physical outcome to the func- 
tional form the harmonic potential, which should provide 
an excellent approximation to any realistic physisorption 
potential, especially one strong enough to allow only rel- 
atively small quantum excursions in the transverse direc- 
tion. 
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FIG. 6: (Color online). Integrated pair correlation function 
computed at T=1.0 for a "He film of coverage 0=0.0719 A'^ 
(solid line), 0.0788 A"^ (dashed line) and 0.0918 A"^ (dot- 
ted line), using the harmonic approximation (Eq. [3]) for the 
substrate-adatom interaction. The two sets of data refer to 
two values of the parameter a, yielding the root-mean-square 
excursions of free atoms in the transverse direction. 



In order to obtain an independent verification of the 
predictions of Ref. (SJ, we have carried out finite tempera- 
ture calculations for the harmonic model ([3]) of confining 
potential, using the same values of the parameter a in 
Eq. ([3]), as well as of coverage 9, as in Ref. |^. Figure [S] 
displays integrated pair correlation functions for two dif- 
ferent values of the parameter a of Eq. [3] and the three 
different coverages. Results shown pertain to simulations 
at a temperature T=1.0 K, but calculations performed at 
lower temperature (down to 0.1 K) yield indistinguish- 
able results, within statistical uncertainties. 

Two different regimes can be straightforwardly identi- 
fied, which correspond to the same already observed us- 
ing the "3-9" potential. Specifically, for the greater value 
of a, corresponding to looser transverse confinement, the 
helium films do not display a 2D character. Indeed, as 
shown by the curves in the upper part of Figure [6l the 
pair correlation functions all remain finite in the region of 
the hard core of the interatomic potential, signalling par- 
ticle layering. The superfiuid fraction in all three cases is 
essentially 100% at T < 0.5 K, i.e., the system features no 
broken translational invariancei^ No hint of crystalline 
structure can be seen. We conclude that the system in 
this case forms a multi-layered liquid, as condensed he- 
lium invariably does on weak substrates. In other words, 
the physical behavior observed here closely reproduces 
that seen on weak substrates, such as Mg. All of this is 
in agreement with the findings of Ref. |8| for this value 
of (7, for which authors find that the ground state is a 
structureless liquid. 

The physics is clearly very different for the case of 
tighter confinement, represented by the lower value of 
a, for which the behavior of the system is essentially 2D, 
atomic motion in the transverse direction being strongly 
suppressed. Indeed, the integrated pair correlation func- 



tion closely resembles that of the purely 2D system, as we 
have verified by simulation. For the two higher coverages, 
the g{r) displays regular, robust oscillations characteris- 
tic of the crystalline phase; concomitantly, the value of 
the superfiuid density is zero, as permutations of indis- 
tinguishable ''He atoms are exceedingly infrequent, down 
to the lowest temperature explored here. Ref. M reports 
minuscule (~ 10"'*) but finite values of the superfiuid 
fraction at T=0, for the same values of a and 9 con- 
sidered here. If those numbers could be regarded as re- 
liable bulk estimates, they would point to a supersolid 
phase for the two larger coverages, with a tiny super- 
fluid response.^* However, it is far from clear whether 
the computational methodology that yielded these esti- 
mates (Diffusion Monte Carlo) truly affords that kind of 
precision. The robustness of these quoted values of the 
superfiuid fraction against various sources of systematic 
error (finite size of the system, population size and trial 
wave function bias, time step error), whose magnitude 
could easily dwarf that of the reported superfiuid signal, 
is not discussed at all in Ref. [^ 

For the lowest coverage, namely 0.0719 A~^, the sys- 
tem is a quasi-2D superfiuid, with no diagonal long-range 
order; the superfiuid density in this case is finite, and sat- 
urates to unity in the low temperature limit. There is no 
evidence that the system may crystallize at sufficiently 
low T, for unlike the case of weak substrates explored 
above, the liquid character here is not associated to the 
formation of multiple layers. Our numerical prediction 
for the lowest coverage case, therefore, is in disagreement 
with that of Ref. |^, wherein a crystalline ground state at 
all coverages is predicted^^ for cr=0.5623 A. 

A few remarks are in order. First, the most recent 
study*^ of the low temperature phase diagram of 2D ''He 
(based on a slightly more attractive version of the Aziz 
potential and smaller system sizes than those utilized 
here) , places this value of coverage inside the very narrow 
coexistence region of liquid and solid phases; this clearly 
renders problematic the unambiguous identification of a 
single phase, especially one displaying two kinds of or- 
der. One would expect, however, that in the presence of 
competing liquid and crystalline phases of very similar 
energy in 2D, transverse fluctuations should strengthen 
the disordered phase. Moreover, the prediction of Ref. 
13 is based exclusively on the resolution of an expectedly 
tiny energy difference (values are not given in Ref. |8|) . It 
is worth reminding that finite temperature methods are 
unbiased, i.e., they require no a priori physical assump- 
tion. On the other hand, the values of physical quantities 
computed with ground state projection methods such as 
the one used in Ref. 8, including the energy, are inher- 
ently biased by the use of initial trial wave functions. 



IV. CONCLUSIONS 

We have carried out extensive numerical studies of ''He 
films adsorbed on weak substrates, searching for a pos- 



sible quasi-2D supersolid phase, i.e., one displaying si- 
multaneously diagonal and off-diagonal long range order. 
We have restricted our search to the case in which density 
long-range order occurs through spontaneous breaking of 
translational invariance, i.e., it is not caused by an exter- 
nal pinning potential. 

The hypothesis tested here is that quantum fluctua- 
tions of atoms in the direction transverse to the plane 
may enhance the tendency of identical particles to ex- 
change, possibly leading to the global superfluid phase 
coherence, while preserving density long-range order. 
The results shown here, obtained for two different models 
of substrate adsorption, constitute strong evidence that 
such mechanism does not take place; only non-superfluid 
crystalline phases or superfluid liquid ones are observed, 



depending on the weakness of the substrate and on the 
temperature. It seems, therefore, that this physical ap- 
proach, namely effective reduction of dimensionality, is 
not likely to lead to the observation of a supersolid phase. 
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